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We present the alternative way of derivation of the Wandzura-Wilczek re- 
lations between the kinematical twist-3 and twist-2 functions, parameterizing 
hadronic matrix element in two-photon processes ^*tt — > jtt and 7*7 — > vrvr. 



^ ■ The new equations, providing the independence of the physical cross-sections 

on the choice of the light-cone direction, are suggested and explored. The 



amplitude of 7*7 — > vrvr up to genuine twist-3 accuracy is found. 



Introduction. The investigation of the two-photon amphtudes of Deeply Virtual 
Compton Scattering (DVCS) or virtual photon production of meson pairs, up to 
1/Q corrections (twist-3 terms) provides one of the interesting lines of inquiry. In- 
deed, these corrections ensure the electromagnetic gauge invariance of the DVCS 
amplitude with the sizeable transverse components of the momentum transfer. This 
problem was discussed in the pioneering paper [|l| and the phenomenological solu- 
tion was presented in [|]. The field-theoretical treatment was suggested in P, H and 
developed in 0, §|, From the other side, the study of these corrections provides 
the possibility to consider the delicate effects of factorization [Rl 0. 
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It was shown that, if the dynamical (genuine) high twist terms due 

to the physical transverse components of gluonic fields are neglected, the param- 
eterizing functions (the skewed, or generalized, parton distribution) of kinematical 
twist-3 can be related to the parameterizing functions of twist-2 by means of the 
Wandzura-Wilczek (WW) relations. Their general structure was investigated, and 
it was demonstrated that within the WW approximation the skewed twist-3 distri- 
bution functions have the discontinuities at x = ±.^, which fortunately do not lead 
to the observable consequences at the level of twist-3 accuracy. However, the general 
relevance of these effects (and, in particular, their possible appearance in the gen- 
uine twist case, for both DVCS and meson pair production channels) is still obscure. 
To some extent, this may be due to the lack of the physical interpretation of the 
WW approximation, essentially based on the involved treatment of the light-cone 
operators with total derivatives in the coordinate representation. 

In this note we present the alternative derivation of the WW relations for two- 
photon processes based on the Ellis-Furmanski-Petronzio (EFP) |^ factorization 
scheme. In the framework of EFP scheme, the separation of the kinematical and 
dynamical twist-3 effects is rather physical and amounts to the taking into account 
the terms linear in the intrinsic transverse momentum of partons in hadrons ^. As 
a result, we present the decomposition of parameterizing functions of twist-3 to the 
WW-type contribution and contribution arising from the three-particle correlators. 
We apply this method for the amplitude 7*7 nn. As a result, it is calculated 
keeping both kinematical and genuine twist-3 corrections. We found the analogs of 
the mentioned singularities and show that the physical amplitude at twist-3 level is 
not affected by them. 

DVCS amplitude of pions in WW approximation. As in [Q, we adhere the EFP 
factorization scheme. Neglecting the dynamical twist-3 contribution, which amounts 
to the keeping of the operators with transverse derivatives alone (c.f. IjTUI)) being 
just the Fourier transform of the mentioned intrinsic transverse momentum, we are 
able to write the DVCS amplitude for the pion target in the following form: 

T^l^+Ti.'J = I rfxtr|E^,(xP)r(x)} + 

J dxidx2tr^E^p^{xiP,X2P)ujpp>Tp>{xi,X2)^ (1) 
where Upp' = 5ppi — Up/Pp, and 

Ta^ix) = -J dXe''^^+^^\p'\M^n)MO)\p), 
^This method in the case of DIS was developed by one of the authors a few years ago 
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r^^(a;i,a:2) = - J d\id\2 e'^'^'^^^^'^'^"'^^' 

(pVa(Ain)z 9^' {X2n)M0)\p), (2) 

where p' = P + A/2, p = P - A/2, A = q - q' = -2^P + Aj- and p{p'), q{q') 
are the initial (final) momenta of pion and photon. The T^"^ and T^^^ denotes the 
contributions from the twist-2 and twist-3 amplitudes (c.f. 0]). 
Parameterization of hadronic matrix elements. In terms of the light-cone basis, 
which consists of the vectors P, n (P^ = = and P ■ n = 1) and Ay, the 
parametrization of the vector and axial hadronic matrix elements (other structures 
do not contribute for massless quarks) can be introduced, as in Namely 0, 

{p'\miM^n)\p) = |rfxe-^(-+«)^|//i(x)P^ + i/3(x)Aj|, (3) 
^ 1 

{p'\mi, < ^(An)b) = / (ia;e-^(^+«)^i/f(a;)P^Aj, (4) 

1 

{p'm^h,l,^{\n)\p) = j dxe~'^''+^^^tHA{x)e^^Tp^, (5) 

-1 

^ 1 

(p'|^(0)757m ^{Xn)\p) = J dxe-"^^^^^\Hl{x)P,ep^Tp^, (6) 

-1 

for the quark correlators, and , 

(p'lV'(0)7M^?<(A2ri)^(Ain)|p) = 
1 

j rfxirfx2e-^(^i+«)^i-*(^2-^^)^2P(xi,X2)P^Aj, (7) 

-1 

(p'l^(0)757A.^7<(A2n)V^(Ain)b) = 
1 

j dXirfX2e-*(^l+«)^l-*(^^-^^)"^2D(Xi, X2)P^£pATPn, (8) 

-1 

for the quark-gluon correlators. Here we see the contributions due to the quark 
intrinsic transverse momentum and transverse component of the gluon field, respec- 
tively. In the WW approximation one is keeping only (§, H), neglecting (|^, ||), and 
in this sense it is already present in |^ . The remaining ingredient is to express these 
^here EpA^Pn = epap.yA'^°'P^n'< . 
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functions in terms of the leading twist one, which can be done using the specific 
version of the rotational invariance. 

The n-independence of DVCS amplitude. Using the parameterization (^)-(^ we 
write the expression for (|l]) in the form 

T^l^+Tj^'J = J dxtr!^E,,{x)P^H,{x,0 + 
J c?xtr|E^,(x)A^|i73(a;,0 + 

J dxtTS^E^^,{x)Yl5^'i'HAix,^)epATPn + 

J dxidx2tTl^Ef,p^{xi, X2)P^Hl {x2, O^i^i - ^2)^1 + 

j dXidX2tT^E^py{xi, X2)Pl5^iH'^{x2, O^i^l - X2)(^pATpn- (9) 

The gauge invariance requires that the physical amplitudes are independent on the 
gauge-fixing parameter which, in our case, is the light-cone vector n. As it is also 
used to define the light-cone direction, the same is required by Lorentz invariance, 
which should be implemented keeping the constraints (ra^ = 0,P ■ n = 1). We 
are therefore dealing here with a particular subgroup of the Lorentz group. The 
relevant symmetry may be called generalized rotational invariance, to distingush 
with the standard rotational invariance in a particular rest frame, and to stress its 
specific form in the case of generalized parton distribution, to be discussed below. 
One should recall that the Burkhardt-Cottingham sum rule, known to be intimately 
related to rotational invariance, follows from the WW approximation for polaritzed 
DIS, and may be derived from n independence [p!2| , |10|. Note that it is rather Lorentz 
than gauge invariance which is crucial here, as one may consider the covariant gauge, 
where role of as a gauge-fixing vector is absent Formally, this invariance can be 
written as 

Tp4n + Sn] = Tp4n], (10) 

or, equivalent ly, 

;{r(".^+r(t)} = 0, (11) 



d 



dnp 

Note the essential difference with the (unpolarized) DIS case f^, where dependence 
on n was coming only through the explicit dependence of the appearing tensor 

•^The new terms, proportional to A • n are absorbed to the standard P-exponent of the leading 
twist operators in this case 
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structures. Here, it comes also through the ^ dependence of the matrix elements, so 
that 



drijs drip 2 dS, 

The second term does not appear also in the polarized DIS case |12[ (where the 
role of A is played by the polarization vector S*), as the dependence of parton 
distributions and correlations on {S ■ n), being the analog of ^, is absent due to the 
general property of the linearity of density matrix. It is this extra term, specific to 
the case of generalized parton distributions, which is one of the reasons to use the 
term 'generalized' for this particular version of rotational invariance. 

Note that the relation (p!TD is in some sense analogous to the renormalization 
group equation governing the dependence of hadronic matrix elements on the fac- 
torization scale. Indeed, the factorization scale /i is the unphysical auxiliary param- 
eter required to write the physical amplitude in the factorized form. Analogously, 
n is the auxiliary vector, required to specify the infinite- momentum frame. The 
dependence of the observable quantities on both these parameters, /i and n, must 
vanish. 

Using the Ward identity, the computation of the derivative in ([Tl|) leads to the 
following equation: 



y"rfxtr|E^^(x)p|A^iJ3(x,0 + j dxii^E^^{x)Yl^iHA{x,i)ep/:^Tpp - 
E^u[xi) - E^^{x2) ]P> + 

/J Xi — X2 

f ^ TP f ^\ \ 'iHl{x2,0S{xi -X2) 
Ef,u{xi) - E^^{x2) 7p75 \ epATp/3 = 0. 

/ J Xi — X2 

(13) 

Note that we actually did not use the specific form of the hard scattering amplitude, 
so this relation should be valid for any hard process. Requiring the validity of this 
equation for arbitrary E'^y, we can easily derive the n-independence condition in the 
more simple form: 

dH^{x,0 _ ldH,{x,0 




dx 2 
dHl{x,i) 



+ H3{x,0, 



dx 



HAix,0- (14) 
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QCD equations of motion for matrix elements. Following 0], we write the QCD 
equations of motion 

1 

dy {B^^\x,y) - D^'\x,y) - 5{x - y)Hl{y) 



-^Hs{x)-^Hi{x)-xHa{x) 



dy {B^'\x,y) + 5{x - y)Hfiy) - D^^\x,y)^ 
xHs{x)+^HAix), (15) 



-1 



where symmetrical and anti-symmetrical functions are defined as, 

y) = ^ {B{x, y) ± B{y, x)) . (16) 

Within the WW approximation they can been rewritten as 

Hl{x,0 = ^H,{x,0 + \H^{x,i)+xHA{x,i), 

iff (x,0 = xH^{x,i)+iHA{x,i). (17) 

Note that these equations of course contain the so-called total derivatives, crucial 
for the derivation of WW relations P, ||, 0. 

The WW relations for DVCS process. Calculating the derivative of ( |17|) with respect 
to X and using (|1^), we obtain the system of differential equations: 

dH+{x,0 _ l d,H,{x,0 

dx 2 X + ^ ^ 

dH4x,0 _ l d+H,{x,0 

dx 2 X — ^ ' 

where we introduced the notations: 



(18) 



8 H(r ..dH^ dH^ 

Hiix,0=Hsix,0±HAix,0- (19) 
Taking into account the boundary conditions 

H±{l,O = H±{-l,O = 0, (20) 
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solution of the system (|T8[) can be easily found: 



a; ^ _1 ^ ^ 

(21) 

The theta-functions appearance is due to the fact, that we are using the different 
boundary conditions below and above the singular points x = Expressing the 
if;^!^— functions by making use of (|19D, we have: 

= -le(x > -0 J - ie(x > 1 + 

ie(. < -0 / iy?^^ + ie(. < {) / d,MiM) , 

4 y + i 4 y~i 

4 J y+^ 4 J y-^ 

4 ^ 1/ + ^ 4 J y-^ 

(22) 



-1 



Note that a given solution completely coincide with the expressions presented in ref. 
1^, although our representation for looks a bit more simple. 

Using these results, we can now decompose the i/3(A)— functions to the WW type 
part and genuine part arising from the three-particles correlators. Using (15) and 
(Pip, we obtain the following 

Hs{x,0 = g3^^(x,0 + l^dy^Bix,y,0-Dix,y,0) + 

^^^^ / dy (s(x, y, -0 - Dix, y, -O) (23) 

and 

1 

Ha{x,0 = g]r^(x, - \ J dy (six, y, ~ D{x, y, o) + 



-1 



2(x + 



1 

J dy(^B{x,y,-0-D{x,y,-0). (24) 

-1 



Here we take into account the properties of B— and D— functions following from 
the T-invariance, i.e. 

B{x, y, = B{y, x, -0, D{x, y, = -D{y, x, -0- (25) 

We can see that the functions and H (both in WW approximation and 
beyond) possess the singularity at points — ^ and respectively. Consequently, as 
it was mentioned in |Q] and confirmed in the functions and have finite 
jumps in these points. One can add here, that if the function Hi{^,^) is finite, this 
singularity is in fact logarithmical, and the mentioned jump means nothing more 
than the finite difference between two logarithmically growing values H^^a{x ± e). 

From the other side, in the DVCS twist-3 amplitude we are dealing with the 
following combinations of if+ and iJ_ : 

l^iP, + (h,{x, + Iha{x, , (26) 

which reduce to the structure at x = ^ and if .—structure at a; = — ^, and 

(^^P. + Q?j (^H,{x, - Iha{x, , (27) 

which reduces to the if„— structure at a; = ^ and structure at x = — ^. 

In (p6[) and ( ^7]) the indices /i and u belong the virtual and real photons, respec- 



tively. Therefore the DVCS amplitude with the longitudinal polarization of virtual 
photon (^) is free from the divergences. At the same time, in the case of the trans- 
verse polarized virtual photon (|27|) the contribution from twist-3 amplitude to the 
observables is suppressed as order after the contraction with the polarization 
vector 0, ||. 

77* —>■ TTTT process. We are now passing to the consideration of — ^ 
7r(p')7r(p) process, adopting the new notations for P = p + p' = q + q' and A = 
p' — p = — + At- Acting as previously, we can derive the relation which is 
the result of independence on the choice of the gauge vector n. All the stages of 
derivation are the same up to the parameterization of the relevant matrix elements, 
known as Generalized Distribution Amplitudes (GDA) |]15|. We are suggesting here 
the natural set of twist-3 GDA consisting of the following ingredients: 

(p,pV(0)7M^(An)|0) = I dxe'^4H,{x;^,s)P, + Hsix;^,s)AA, (28) 
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1 



(p,pV(0)757/.^(An)|0) = J dxe'^^iHA{x;^,s)e^^Tp^, (29) 



^ 1 

(P,PV(0)7^^ i^{Xn)\0) = J dxe''^H^{x; s)P,A^^ (30) 



} 

{p,p'\mW dj V(An)|0) = J dxe'^^iH^{x;C,s)P^e,^Tp„ (31) 



1 

J dxidx2e''''^'+'^'''-'''^^'B{xi,X2)P^A^p, (32) 



{p,p'\mi5lM^M2nmXin)\0) = 

1 

y ciXiCiX2e'^i^i+*(^^-^i)^^i£'(Xi,X2)P/.£pATPn- (33) 


In terms of these parameterizations, the n-independence condition for a given case 
reads as 

'-^^ = HAM), (34) 

while the QCD equations of motion take the following form 
1 

/ dy {m{x, y) + 5{x - y)Hl{y) - B^^^x, y)) = 



^Hs{x,^)+(x-^^Ha{x,0: 
1 

/ dy{m'\x,y) + 5{x-y)Hfiy)-mix,y)) = 



I^x-^^Hs{x,^) + ^Ha{x,0: (35) 



so that within the WW approximation we have: 



1\ . 1 



H^ix, = - 2) Hsix, + ^H^ix, 0- (36) 



Differentiating (^) with respect to x and taking into account ([3^), we can reduce 
the system of equations (|3^ and (^) to the following form ^ 



) 

(37) 



dx X dC, ^ 

dH-{x,0 1 dHi{x,0 



dx {x — 1) 

Note that we cannot now always implement zero boundary conditions for both 
functions due to the singularities in the r.h.s. The only possible zero boundary 
conditions are: 

^+(1,0 = 0, 

^-(0,0 = 0. (38) 
It can be easily found that the solution takes the form ^ 

^r'^(x,0 = -/rfy^^^, (39) 

y 

X 

H^^{x,0=Jdy^-^^^, (40) 

or, for H^^— functions: 

H'n-, f ) = i / dy?i^ - \ j , (41) 



HrM--,jiv^j^--,jiv^^^ (42) 







These expressions are analogous to WW relations for the transverse polarized p- 
meson light-cone distribution derived by Ball and Braun [|l^. This is not surpris- 
ing, as the momentum just plays the role of the polarization vector ex- In 
particular, the structures if+ and correspond to the distributions with a defi- 



nite (anti)quark spin project ions |T^, whose zero endpoint values result from their 
conformal properties. 

^ H±— functions are defined in complete analogy with H±— functions, see ( pj| ) 
^here = d/d^ 
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We can now easily obtain the decomposition the if3(A)— functions on the WW 
type contribution and the three-particles contribution, as it was done above for the 
DVCS case. We have 

Hsix,0 = H^'^ix, - \ . I dy (^Bix, y, - D{x, y, o) + 

^ldy(^Bix,y,0 + Dix,y,0) (43) 



and 



HAix,0 = g7^(x, + _^J dy {b{x, y, i) - D{x, y, o) + 

^ ' 

^jdy{B{x,y,i) + b{x,y,i)y (44) 



In ( |43D and (|4^), we used the properties of C-invariance which are the following 

B{x, y, = B{y, x, 0, D{x, y, = D{y, x, 0- (45) 

We here stress again that in both cases, i.e. in the genuine expressions (^), 
( |I^ ) and WW approximation (|41|), (^2]), the functions if+ and if_ have, generally 
speaking, singularities at the points x = and a; = 1, respectively. However, if 
one assumes zero boundary condition for function Hi, required by the factorization 
at the leading twist level, these singularities are transformed to the non-zero finite 
limits at these points. They are the counterparts of the x = ±^ singularities in the 
DVCS case. Anyway, these finite values would be sufficient to spoil factorization 
after the contraction with the short- distance subprocess cross-section. As in the 
case DVCS process, this effect can not affect the amplitude contracted with the 
polarization vector of the real photon. 

Indeed, we can derive the expression for the 77* —>■ tttt amplitude using the 
developed approach |Q. It reads 

1 



(46) 



where 



Vf = -q'gU'^x-l)Hi{x,0 + 
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AX({2x-l)Hsix,0-HAix,0 

Note that, contrary to DVCS amplitude case, the contribution of Hi is gauge invari- 
ant by itself. This is due to the fact, that parameterization of (anti) quark momenta 
in terms of the vector P makes both of them coUinear and independent on A^, so 
that one is dealing with the physical on-shell amplitude, which is gauge invariant. 
At the same time, only the total contributions of new quark and quark-gluon matrix 
elements is gauge invariant provided the equations of motion are taken into account. 
The contribution for the longitudinally polarized photon arises from the term 

(P + q'), ((2x - l)^3(x, + Ha{x, o) , (47) 

which is described by the structure at x = and by structure at x = 1. 
Let us recall that at these points H_ and functions are finite. The respective 
contribution to the amplitude reads as 

TS)''*=-\ I ^r^Al{P + q'),[{2x-l)Hs{x,0+HA{x,0)- (48) 
q ^ x{ X) 

Using (|1|) and (|4|), we substitute the H^^- and -functions in Eq. (|8D 

which thereafter can be rewritten in the form 

r«-^-^^<^aJ..(4M),£iM)). ,«) 



where 

Mx,0 = f dy^ii^, M,,o = j dy^ii^. (50) 
After simple manipulations, we derive the following expression 

r»)-- = ^^<^±«%/d,a,ff,(,.{)(i5^^ - (51) 



which coincides with the expression in At the same time, the genuine twist-3 
contribution is new. 

Conclusions. In this paper we have presented the new derivation of the WW rela- 
tions between the parameterizing twist-2 and kinematical twist-3 functions. This 
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simple approach is based on the use of the equations of motion and on the require- 
ment of independence of the amphtudes from the gauge-fixing hght-cone vector. 

We have found the exphcit expression for the genuine twist-3 functions for DVCS 
and 77* —>■ 2tt processes. We also have computed the amplitude of 7*7 tttt 
including complete twist-3 corrections. 

The counterparts of the singularities appearing in the DVCS case are provided 
by the finite values of generalized distribution amplitudes in the endpoints, which 
are related to their conformal properties. Therefore, one may hope that appearance 
of the singularities, as well as their cancellation in the physical amplitudes, have 
some general reasons. 

We have also shown that the contributions from the three-particles correlators 
posses the same singularities as the WW parts which also do not contribute to the 
physical amphtudes at the twist-3 level. 
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